A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 

02.01.1; 02.19.1; 02.20.1; 09.03.2; 09.04.1; 13.07.1 



ASTRONOMY 

AND 
ASTROPHYSICS 



Proton and electron acceleration through magnetic 
turbulence in relativistic outflows 

R. Schlickeiser 1 and C. D. Dernier 2 

1 Institut fur Theoretische Physik, Lehrstuhl IV: Weltraum- und Astrophysik, Ruhr-Universitat Bochum, D-44780 Bochum, 
Germany; r.schlickeiser@tp4.ruhr-uni-bochum.de 

2 E. O. Hulburt Center for Space Research, Code 7653, Naval Research Laboratory, Washington, DC 20375-5352 USA; der- 
mer@osse.nrl.navy.mil 



o : 
o , 
o ■ 



o , 
\o ■ 
vn 

o : 
o , 
o ■ 

On; 

6 ' 
a: 



Received 16 February 2000 / Accepted 15 May 2000 



Abstract. Low frequency electromagnetic turbulence is 
generated in relativistically outflowing plasma that sweeps 
up particles from the surrounding environment. Electrons 
are energized by stochastic gyroresonant acceleration with 
the turbulence produced by the isotropization of captured 
protons and charged dust. Protons are accelerated stochas- 
tically by the dust induced turbulence. Analytical solu- 
tions for the proton and electron energy distributions are 
obtained and used to calculate broadband synchrotron 
emission. The solutions are compared with generic spec- 
tral behavior of blazars and gamma-ray bursts. Dust cap- 
tured by a blast wave can generate turbulence that could 
accelerate protons to very high energies. 

Key words: acceleration of particles - shock waves - tur- 
bulence - cosmic rays - dust - gamma rays: bursts 



1. Introduction 

Measurements of radio and optical polarization of blazar 
jet sources indicate that nonthermal lcpton synchrotron 
emission is the dominant radiation mechanism at these 



frequencies (Begelman et al. |1984|) . The radio through 
X-ray afterglow and the promp t soft-gamma ray emis- 
sion observed (C osta e t al. 1997 , v an Pa radijs et al. 1997 , 
Djorgovski et al. 1997 , Frail et al. 199 j ) from gamma-ray 
burst s (GRBs) is attribute d to t he same process (Tavani 
1991 , Paczyhski & Rhoads |1993| , Meszaros & Rees |1993| ). 



In the standard blast- wave p hysics that has been devel- 
oped (Meszaros & Rees |1997[ Victri |L997|, Waxman [L997D 



to explain the long wavelength afterglow emission from 
GRBs, particles captured from the external medium en- 
ergize the blast wave and cause blast wave deceleration. 
When a blast wave sweeps up material from the surround- 
ing medium, the free energy of the captured particles ini- 
tially resides in the more massive protons, nucleons, or 



charged grains and dust, if the latter can survive capture 
by a relativistic wind, as assumed here. The less massive 
nonthermal electrons and positrons emit, however, most 
of the radiant energy. In order that high radiative effi- 
ciencies are possible, it is therefore required that there is 
efficient transfer of energy from the swept-up massive par- 
ticles to electrons, or to very high-energy protons emitting 
synchrotron radiation or photomeson secondaries. 

Major uncertainties in modeling relativistic outflows 
involve magnetic field generation and the mechanism for 
dissipating and transferring the free energy of protons, 
ions and dust particles to lighter particles (Meszaros, Rees 
& Papathanassiou 1993, Chiang & Dermer 1999). Non- 
thermal electrons and protons may both undergo first- 
order Fermi shock acceleration in relativistic blast waves 
to form power-law distributions reaching to very high en- 
ergies, and the origin of ultra-high energy cosmic rays 
(UHECRs) has been attributed to shock acceleration of 
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protons in the GRB blast waves (Vietri 1995, Waxman 
|ll)95 ). This mechanism has recently been called into ques- 
tion (Gallant & Achterberg 1999] ) . Following the first shock 
crossing by a particle, subsequent cycles do not permit 
large gains of energy because the particle is captured by 
the shock before it has been scattered through a large 
angle. A stochastic mechanism for particle energization 
in relativistic shocks avoids these difficulties, but requires 
strong ma gneti c fields to retain UHECRs within the shell 
(Waxman 1995 ). However, the detailed physics of stochas- 
tic gyroresonant acceleration was not considered. 

In this paper, we show that the magnetic turbulence 
generated by the capture of heavy charged particles can 
accelerate lighter charged particles to high energies through 
stochastic gyroresonant acceleration. This could account 
for the hard spectra observed from GRBs and blazars dur- 
ing their flaring states. The acceleration of particles to 
high energies could produce UHECRs through gyroreso- 
nant acceleration in relativistic outflows. The question of 
the survivability of dust against sublimation by the intense 
radiation fields of the GRB is addressed in the Appendix. 
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2. Particle energization in relativistic outflows 

2.1. The relativistic pick-up model 

A two-stream instability is formed when a blast wave with 
an entrained magnetic field encounters a medium of den- 
sity n *. Ac cording to the recent work of Po hl & Schlick- 
eiser ( 2000 hereafter referred to as PS 2000), this causes 
the captured protons and electrons to rapidly isotropize in 
the blast wave plasma. The isotropization p rocess operates 
on the comoving time scale (eq. (74) of PS 2000| ) t iso (s) = 

1 /2 

0.7n b s /Tsoon*, where T = 300r3oo is the bulk Lorentz 
factor of the blast wave and n&,8 denotes the comoving 
blastwave density in units of 10 8 protons cm -3 . The blast 
wave particle density in the comoving frame is given by 
7ib(cm~ 3 ) ~ E/(Y m p c 2 Airr 2 Ar)~ 10 8 £ , 54/r 3 6 , and is im- 
plied by depositing an initial energy 10 £54 ergs/(47r sr) 
in a blastwave shell with comoving radius Ar ~ r/^o, 
where Tq is the initial blastwave Lorentz factor. Alfvenic 
magnetic turbulence is generated by the incoming protons 



and electrons with an energy density (eq. (53) of PS 2000) 



AUa = f3Am p c' 2 n*T\T 2 - 1, where (3ac is the Alfven 
speed. The total energy density of the swept-up mate- 
rial in the comoving frame is AUi = m p c 2 n^ 
(Blandford & McKee 1976). Thus the fraction of incoming 



1 



energy converted to transverse plasma wave turbulence in 
the process of isotropizing the particles is ~ (3a ■ 

P Sj2000 calculated the emitted radiation in the thick- 
target limit when the captured, isotropized protons un- 
dergo nuclear inelastic collisions with thermal particles in 
the blastwave. Although this process involves the fraction 
~ (1 — Pa) of the incoming energy retained in the pro- 
tons, it operates on the comoving pion production time 
scale t pp (s) ~ 1.4 x 10 7 /n(,,8, where n^g denotes the co- 
moving blastwave density in units of 10 8 protons cnr 3 . 
If the acceleration rate is sufficiently rapid, pion produc- 
tion could produce low-level prompt and extended emis- 
sion in GRBs or blazars if remains at the estimated 
level. Moreover, luminous radiation from pion-decay elec- 
trons and positrons is only produced if t pp is much shorter 
than the diffusive proton escape time scale from the blast 



wave, given by t e 



{Arf/n = r 2 /(T ct iso ) 2 , noting 



that the spatial diffusion coefficient k = c t iso . This re- 



quires that n 



1/2 

6,8 
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» 15P 



300 



which holds well in jets 



but is only marginally satisfied in GRB blast waves. 



2. 2. The role of captured charged dust 



Here we extend the work of PS200C by considering the 
process of channeling the plasma turbulence energy into 
the energy of the swept-up primary electrons through gy- 
roresonant interactions. As demonstrated in that paper, 
incoming charged particles with mass rrij, charge Qje, and 
Lorentz factor T generate Alfvenic turbulence with par- 
allel wavenumbers Ifel > -R7 1 , where the Larmor radius 



Rj — rrijC 2 \/T 2 — l/(Qj\e\Bo) and Bq is the entrained 
magnetic field in the frame of the blastwave. For simplicity, 
we assume that the magnetic field is parallel to the blast 



wave velocity, and generalize the calculation of P52000 
by including the effects of negatively charged dust. Ac- 



cording to eqs. (48)-(50) of PS200C, we obtain the power 
spectrum of backward (— ) and forward (+) propagat- 
ing Alfven waves given by I—(k) = |Z(A;)| and I+{k) — 
I 2 (k)/\Z(k)\, where Io(k) is the initial turbulence spec- 
trum before isotropization, and 

(k) _ 1 Up4 n p B 2 



E 



p,z ""p 

c n b 
Qj rij 



n 



[1 - {Rjk)- 1 ] H[ \k\-BZ 



(1) 



v 



The term cjp j(s _1 ) = 1.3 x 10 3 y / n& is the proton plasma 
frequency and H[x] denotes the Heaviside step function. 
This expression assumes Io(k) <C |Z(fc)| and makes use of 
the charge-neutrality condition QdV-d + n e = n p = Tn*, 
where we assume that all dust particles have the same 
charge. A number of processes influence the charging of 
dust grains in a proton-electron-plasma. In the absence 
of any other process due to the higher electron mobility 
an initially uncharged dust particle will absorb more elec- 
trons than protons and thus attain a negative potential of 
order — 2.5fcsT e /e (Spitzer 1968). However, a number of 
competing processes like neutral-atom collisions and ion 
impacts on the grain surface lead to electron loss from 
the dust grains. Also, a significant ultraviolet photon flux 
can charge the grain by the photoelectric effect to a pos- 
itive potential. All these processes, together with poorly 
known grain properties, make an exact determination of 
the grain charge impossible (see e.g. McKee et al. (1977) 
for the variety of grain properties resulting from model 
calculations under various assumptions). Here we follow 
the argumentation of Ellison, Drury & Meyer (1997). If 
the grain potential is 4>, then the charge on a spherical 
(of radius a) grain is of order q = eQd — 47reoa</>, so that 
Q d ~ 7OO(a/(lO" 7 m))(0/lOV). The number of atoms in 
the grain will be of order [a/(10 _10 m)] 3 , so that with the 
mean atomic weight /i of the grain atoms the entire grain 
atomic weight is Aq = /j[a/(10~ 10 m)] 3 . This implies a 
dust mass of = fim p [a/ (10~ 10 m)} 3 = 2 • 10 10 m p a 3 _ 7 
if we adopt fi — 20 for silicate grains. Combined with 
the charge estimate this yields for the charge-weighted 
dust/proton mass ratio Yd = fnd/{Qdjnp) — 3 ■ 10 7 a 2 _ 7 . 

It is important to note that the densities of the incom- 
ing protons, electrons and dust particles are much smaller 
than the density of the blast wave plasma nl , so that the 
properties of the low-frequency Alfven waves, carried by 
the blast wave plasma, are not modified by the captured 
particles. Mathematically, this corresponds to the state- 
ment that the real part of the Alfvenic dispersion relation 
is solely determined by the blast wave electrons and pro- 
tons. However, the free energy of beams of incoming pro- 
tons, electrons and dust particles gives rise to a non-zero 
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imaginary part of the Alfvenic dispersion relation lead- 
ing to a positive growth rate of backward moving (in the 
blast wave plasma) Alfven waves. Moreover, for relativis- 
tic beams (Miller 1982 ) the growth rate of longitudinal 
electrostatic waves is much smaller than the growth rate 
of transverse Alfvenic waves, so that the free energy is dis- 
sipated by isotropising the incoming beam (i.e. capture) 
rather than heating the blast wave plasma by plateauing. 

Two comments on the equilibrium wave spectrum (1) 
are appropriate: 

(1) According to the quasilinear relaxation theory the 
equilibrium Alfven wave spectrum given in Eq. (1) is achieved 
formally after infinitely long time t^i — > oo by transferring 
the free energy in the initial particle-beam-distributions 
into plasma waves. Numerical simulations of the electro- 
static beam instability (see e.g. Grognard 1975) indicate 
that this asymptotic equilibrium distribution is established 
after t q i — 100ti SO . 

(2) The wave spectrum in Eq. (1) is calculated from quasi- 
linear wave kinetic equations that only take into account 
the wave growth from the unstable particle beams but ne- 
glect wave-wave interactions. Since the latter scale with 
the total magnetic field fluctuation energy density (SB) 2 , 
this is justified as long as the total wave energy density 
is small to the initial energy density of the beam parti 



10' 



10° r 



o 

CM 



CD 



cles. Because according to Eqs. (53) and (64) of PS |2000 
(SB) 2 /(4ttwo) = Pa « 1 this is indeed fulfilled. 

The turbulence spectrum (1) leads to the quasilinear 
isotropization time scale ti so given above. This time scale 
and tq\ = 100ti so are much shorter than the blastwave 
crossing time i/ c (s) ~ Ar/c ~ 1000ri 6 /r 30 o. Fig. 1 shows 
the resulting turbulence spectrum of backward propagat- 
ing, left-handed circularly polarized Alfven waves for a 
charge- weighted dust/proton density ratio e = Qdrid/n p = 
0.01, and a charge-weighted dust/proton mass ratio Yd 
= Trid/(Qd'm'p) — 10 6 . The relative contributions of dust 
to the turbulence spectrum (1) is determined primarily 
by the parameters e and Yd- Most of the magnetic turbu- 
lence energy density resides at low wave numbers and is 
provided by the charged dust for this value of e. 

The total dust mass Md = rridUd = YdeM p , where the 
total proton mass M p = m p n p . Thus the total mass in dust 
for the calculation in Fig. 1 is ~ 10 4 times the mass in the 
ionized gas, as might be found in the weakly ionized, low 
temperature star-forming environments that harbor GRB 
sources. By contrast, the relative mass in dust to that in 
protons would be much less in the enviroments of active 
galactic nuclei, where the intense UV and X-ray radiation 
field would photo-disintegrate the dust. 

2. 3. Stochastic gyroresonant acceleration of electrons and 
protons 

The primary electrons and protons may increase their en- 
ergy by stochastic gyroresonant acceleration with the tur- 
bulence spectrum given by Eq.(pJ). Since a fraction (3a 
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Fig. 1. Intensity I-(k) of backward propagating Alfven 
waves at wave vector k generated by a two-stream instabil- 
ity when charged dust, protons, and electrons are captured 
by a relativistic blast wave. Parameters: e = Qdnd/n p = 
0.1 and Y d = m d /(Q d m p ) = 10 6 . 

of the incoming free energy is converted into transverse 
plasma wave turbulence, these particles can tap at most 
a fraction Pa of the swept-up energy via this process. Ef- 
ficient energy transfer therefore requires that (3a ^ 0.01. 
The e-folding time scale for stochastic gyroresonant accel- 
eration is t acc ~ (3(^ 2 ti SO: whereas the d iffusive escap e tim e 
t csc ^ R 2 /n = t 2 JU so (Barbosa [l979| , Schlickeiser |1989[ ). 
Because the turbulence spectrum in Eq.(Q) is oc k~ 2 , the 
acceleration, escape, and isotropization time scales are in- 
dependent of particle energy. The comoving deceleration 
time scale t dcc (s)= 1.3 x 10 i (E 54 /n*Tl 00 ) 1 / 3 (Meszaros k 

Rees |1993D . Provided that (3 A > 0.03n^ 8 4 (r§ 00 /£ 5 4< 2 ) 1/3 , 
the acceleration time scale is much shorter than the time 
scale for evolutionary changes of the blast wave due to 
bulk deceleration. In this steady-state calculation 

of the equilibrium particle energy spectrum is justified, 
which results from the balance of gyroresonant accelera- 
tion, diffusive escape, deceleration and radiative losses. 

The steady-state kinetic equation for the phase-space 
density f e of relativistic electrons is (Schlickeiser 1984) 

)fe]-^ = S ( M- 



1 



d_[Z 

7 2 t; 

where S e 



4 C 



,_r 



7 



syn 



F) (2) 



= Vr 2 — ln*c/(47rrr) and A is the area of the 
blast wave that is effective in sweeping up material. 

If no dust is present, then the electrons can only be 
accelerated to 7 < TR p /R e = Tm p /m e . In this range, 
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t a = t acc and tg = i csc . If dust is present, then electrons 
can be accelerated throughout the range Tm p /m e < 7 < 
7cmax = Firid/(rn e Qd), but the acceleration and escape 
times are modified according to the relations t a — t acc / e 
and tE — ricsc- The radiative loss time scale is i syn (s) = 
67rm e c/ '{(J T Bl) = 7.7 x lO 8 ^ 2 . 

Likewise, the steady-state kinetic equation for the phase- 
space density f p of relativistic protons is 
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where S p = ^JT 2 - ln*c/ (AnTr). Eq.(|) holds in the pres- 
ence of dust when 7 < 7 P!max = TR d /R p = Tm d /(m p Q d ) 
= TYd, and assumes that effects of photopion production 
do not limit the acceleration of protons to the highest en- 
ergies. Evidently, the parameter Y d and the bulk Lorentz 
factor r determine the maximum proton Lorentz factor 
7p,max- If no dust is present, then protons are not signif- 
icantly accelerated by gyroresonant processes unless long 
wavelength MHD turbulence is generated through pro- 
cesses not treated here. 

The exact solution to Eq.(||) for electrons accelerated 
by turbulence with a spectrum given by Eq.(Q) is 

f[fi - (3 - o)/2] 



o 



10 h 



, electrons 



10 



iog 10 y 



/ e (r < 7 < 7 0jmax ) = t a s e 

p/i+(a+l)/2,yU-(a+3)/2 



r[i + 2fj] 



2/1 

7* 



exp(-7/7s)x 



M( M . 



r 



1 + 2/i, —)U{n 



1 + 2/x,^) (4) 



(Schlickeiser ]1984| ), where fi = ^/A + (3 — a) 2 /4, A = t a /t E , 
a = t a /tdcc, r(ao) denotes the gamma function, and M{a\, ai 
and U(ai,a2,a,3) are the Kummer functions of the first 
and second kind, respectively. The value 7 S = t syn /t &cc = 



Fig. 2. Steady-state distribution of electrons and protons 
for 7 > T = 300 resulting from stochastic gyroresonant 
acceleration with waves generated in a relativistic blast 
wave for parameters given in text. 



In Fig. 2, we show the resulting proton and electron 
spectra, choosing r 30 o = 1, (3a = 0.03, ri 6 = 0.1, n^g 
1 1 . 1 . e = 0.1, Y d = 10 6 , and n* — 1 cm~ 3 . As can be seen, 
the electron spectrum dN e /d^ oc 7 _1 exp(— 7/11700) is 
very flat below the synchrotron cutoff, whereas the pro- 
ton spectrum displays a spectrum dNp/dj oc 7 -1 - 36 be- 
the maximum Lorentz factor 7 Pimax = 3 x 10 8 . For 



3/2 



is the electron Lorentz factor where the 



370r 3 oon*/n b8 

electron synchrotron loss time scale equals the accelera- 
tion time scale. In the regime 7 > T, an accurate approx- 
imation is obtained by taking the asymptotic expansion 
of the function U for small arguments, resulting in 

/e(T < 7 < 7o,max) = 



2/j. 



r M+(l+a)/2 -(3+a)/2-p, 



exp(-7/7 s ) 



(5) 



The exact solution to Eq.(|J) for protons accelerated 
by turbulence with a spectrum given by Eq.(|l]) is 

fp(T < 7 < 7p, m ax) = 



^acc Sp 

2e/z p 



pp. p + (a p + l)/2 _^,-(a p +3)/2-p. p 



(6) 



where [i p = yj\ p + (3 - a p ) 2 /4, A p = t acc /e 2 t csc , and a p 



these parameters, we find comoving acceleration times of 
~ 330(/3 y i/0.03) _2 s for electrons, and ~ 10 times larger 
for the protons. The choosen value for (3a implies a mag- 
netic field strength of 14 Gauss. Larger values of Pa will 
shorten the acceleration time scales and improve the ac- 
celerated particle efficiency. Requiring that the gyroradii 
of the swept-up dust particles R d = Y d R p in a mag- 
netic field of strength \ABu Gauss is smaller than the 
comoving shell radius 10 14 ri6Lj~ 1 cm limits the charge- 
weighted dust/proton mass ratio to values smaller than 
Y d < 5 ■ 10 6 tiq B14 r]~ Q, and consequently the maximum 
proton Lorentz factor in the comoving blast wave frame is 
7 P ,max = ry d < 5 ■ 10 s r i6 B 14 , IYoq. 

Fig. 3 shows the resulting synchrotron emission spec- 
trum produced by the electron distribution function (5). 
For the parameters used here, the synchrotron spectrum 
is very hard. Particle spectral indices dN/dj oc 7~ p with 
p > 1 can be produced through this process, leading to 
synchrotron spectra F v oc v~ a with a = —1/3, corre- 
sponding to the index of the elementary synchrotron emis- 
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in relativistic blast waves through the first-order Fermi 
mechanism. 
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Fig. 3. Synchrotron emission spectrum in a vF v represen- 
tation for the electron distribution given in Fig. 2. 



sivity formula for monoenergetic electrons. In standard 
shock acceleration theory (Blandford & Eichlcr 1987), par- 
ticles are accelerated with p > 2, leading to values of 
a > 1/2. Harder synchrotron spectra can however be pro- 
duced through nonlinear effects in first-order Fermi ac- 
celeration, or by introducing a low-energy cutoff in the 
electron distribution function. Synchrotron spectra from 
flarin g blazar so urces (Catanese et al. 1997 Cohen et al. 
1997 , Pian et al. 1998| ) and GRBs display very hard spec- 
tra requiring nonthermal lepton spectra with p < 2, con- 
sistent with the stochastic Fermi acceleration mechanism 
proposed here. 

The hardening above the ankle of the cosmic ray spec- 
trum and beyond the ZGK cutoff could be related to a 
hard UHECR component formed throug h sto chastic ac- 
celeration in GRB blast waves ( Waxman 1995 , Takeda et 
1998|) . Because turbulence is generated in the process of 



al 

capturing and isotropizing particles from the external en- 
vironoment, a wave spectrum that can yield very high en- 
ergy electrons and protons through gyroresonant accelera- 
tion must be present in the relativistic blast waves. If dust 
survives capture by the magnetized plasma wind, then 
long wavelength turbulence can be generated that could 
accelerate protons to ultra-high energies on time scales 
shorter than the duration over which the blast wave decel- 
erates to nonrclativistic speeds. Gyroresonant acceleration 
of particles by MHD turbulence provides a mechanism to 
transfer energy from heavier to lighter species, and over- 
comes difficulties associated with accelerating UHECRs 



3. Summary and conclusions 



We have shown that efficient and rapid stochastic gyrores- 
onant acceleration of electrons and protons occurs when 
interstellar protons and dust are captured by a relativis- 
tic outflow. The MHD plasma wave turbulence, created 
within the blast wave plasma by the penetrating pro- 
tons and charged dust particles, is transferred from the 
heavier to the lighter particles through the gyroresonant 
acceleration process. This mechanism is particularly at- 
tractive for the energization of nonthermal electrons in 
GRB sources, as radiation modelling often requires ap- 
proximate energy equipartition between electrons and pro- 
tons (e.g., Katz 



1994 , Beloborodov & Demianski [1995 , 
Smolsky & Usov |gg§ Chiang & Dermer |1999| ). Capture 



by and isotropization within a blast wave, as discussed 



by Pohl & Schlickeiser (200C), provides primary electrons 
with an energy 1836 times smaller than that of the pro- 
tons. As shown here, the subsequent acceleration of elec- 
trons in the proton-induced turbulence leads to approxi- 
mate energy equipartition of these nonthermal electrons. 
Moreover, if in addition charged dust particles penetrate 
the blast wave plasma, we find that further acceleration 
of the particles is possible to an energy determined by the 
radiative loss processes of the electrons and protons. For 
some parameter regimes, this could accelerate particles to 
~ 10 20 eV and account for UHECR production. 

By solving the appropriate kinetic equations, we cal- 
culated the resulting particle energy spectra of the ac- 
celerated protons and electrons, and the observable syn- 
chrotron radiation spectrum of the accelerated electrons. 
The latter peaks at hard X-ray energies and is very flat, 
in agreement with observed photon spectra of GRBs and 
some blazars during their flaring states. 

A crucial issue, particularly for the stochastic acceler- 
ation of protons, is the penetration of dust. While large 
amounts of dust are present in the environments of star- 
forming regions in galaxies, it is also necessary that the 
dust not be destroyed during capture by the blast wave or 
by the intense radiation associated with captured nonther- 
mal particles in the blast wave. We address some aspects 
of this issue in Appendix A. 
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4. Appendix A: Survivability of dust in GRB 
radiation field 

A major uncertainty in this study is whether the dust can 
survive the intense GRB radiation field by sublimation 
before the blast wave arrives and sweeps up the dust par- 
ticle. In the inhomogeneous environment surrounding a 
massive star/GRB progenitor (Dermer & Bottcher 2000| ), 
it is probable that the blast wave will impact a dense dusty 
cloud after passing through a nearly evacuated region. In 
that case, there will be little radiation in advance of the 
blast wave, because the blast wave will only start to be 
energized and radiate after it encounters the cloud. Under 
such circumstances, the survivability of a large fraction of 
the dust to sublimation is assured. 

In the ideal case where a GRB blast wave passes through 
a uniform density medium before intercepting dust par- 
ticles, we outline a calculation of dust sublimation and 
identify parameter sets that are compatible with dust sur- 
vivability. We follow the treatment of Waxman & Draine 



(1999) for dust sublimation. Although reverse shock emis- 



sion could be important for dust sublimation, we neglect 
it here. It has only been detected from GRB 990123 and 
can be strongly suppressed if, for example, the blast wave 
shell is thin and the reverse shock traverses the shell before 



becoming relativistic (Sari et al. 1996 ). 

In this case the optical-UV emission irradiating the 
dust, which is the photon energy range most important 



for dust sublimation (Waxman & Draine 1999), can be 
described by the L v = Lr/, ma x(* / /* / o) 1/ ' 3 portion of the syn- 
chrotron emissivity spectrum where the quantity L u max 



47m*£ 3 m e c 2 crTri?o/9e (Sari et al. 1998), where e is the 



proton charge. Using standard blast-wave physics, B (G) 
where es < 1 is the magnetic- field pa- 
rameter. To give good fits to the prompt phase of GRBs, 
it is necessary that the electrons are in the weakly cooled 
regime (Chiang & Dermer 1999, Dermer et al. 20001 ). For 
a weakly cooled spectrum, i/ Q — v m — Tj^eB /27rm e c, 
where j m = e e T(p e ~ 2)m p /[m e (p e - 1)] and e e >~ 0.1. 
Taking the electron injection index p e i=s 2.2 - 3.0 = 2.5, 
as implied by burst spectroscopy in the prompt and after- 
glow phas e, we find ^ m (Hz) =4.1 x 10 n e 2 r 4 y/nfes (see 



Sari et al. |l998[ Dermer et al. [2000| for more details). The 
luminosity (ergs s -1 ) between the 1 cV = 2.41 x 10 14 Hz 
and 7.5 eV band is 



Li- 



7.5 



1.28 x 10- 5 < 4/ Wee) V Vr 2 / 3 



(7) 



2.0 X 10 46 (^) 4 / 3 (^)^(^)-/3 ( _^)3 ( _L )2 /3 



io- 



10 16 cni v 300 ; 



which also defines our standard parameter set. 

The relationship of T and x to observer time t for an 
adiabatic blast wave decelerating in a uniform surrounding 
mediuip- 



where the deceleration radius Xd — (SEo/SnTlmpC 2 !!,*) 1 ^ 3 
= 2.1 x 10 16 (i?54/r| 00 n2) 1 / 3 cm, the dimensionless time 
t = t/td = ctT^c/xd, and the apparent isotropic explosion 
energy is 10 54 i?54 ergs. 

For the sublimation of dust, we consider the total en- 
ergy change due to (1) heating by the GRB in the 1 - 
7.5 eV band, (2) thermal reradiation by dust, and (3) the 
energy carried away by the sublimed particles. In this ap- 
proximation, the heating rate dT/dt due to dust energiza- 
tion is given through 

_ _ _ L i-7A t)a 2 (t) 

dt 



[4 



3 , , p m_i<£T dE 



dt 



Aitr 2 

4r:a 2 (t)a SB T 4 - 12rrk B T(^)a 2 (t) • ^ . 

m dt 

The rate of shrinkage of dust grain radius a(t) with time 
t, is given by 

§ = -(-) 1/3 -oexp(-6ABT), 
dt p 



(9) 



(10) 



(Waxman 
1 x 10 15 



Draine 



1999), where p/m = 10 23 cm 



F 



) Vq — 

s" 1 , and b/k B = 7 x 10 4 K for refractory grains. 
We look for solutions of Eqs. (pi)- (jl(J) where the dust 
grain is closer than the blast wave to the GRB source 
at the moment the dust becomes completely sublimed. 
There are no interesting solutions for our standard pa- 
rameters; in this case, dust at 10 16 cm becomes sublimed 
when the blast wave has only reached x = 1.6 X 10 15 cm. 
It is easy, however, to find solutions for slightly less ener- 
getic GRBs in more tenuous environments. For example, 
the blast wave will reach dust particles before they sublime 
if the dust resides within x = 9 x 10 15 cm when Eq — 10 52 
ergs and n* = 1 cm -3 . These values of the energy and 
density are not too different than the values found when 
mode ling the afterglow of GRB 980425 (Wijers & Galama 
|1999D . Thus we find that dust will survive destruction by 
sublimation in some, though not all, GRB environments if 
the surrounding medium is uniform. A clumpy surround- 
ing medium improves the chances of dust survival. 
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